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^^ ' Abstract: We further discuss the QCD Reggeon field theory (RFT) as it emerges from the 

JIMWLK/KLWMIJ evolution equation and beyond. We give an explicit expression for the calcula- 



\^ • tion of scattering amplitude in terms of the eigenstates of the RFT Hamiltonian. We point out that 

^"-"^ the spectrum of RFT is doubly degenerate, the degeneracy being related to the spontaneous breaking 



of the Dense-Dilute Duality symmetry of RFT. The degeneracy is between the "almost white" states 
(the Yang sector) which contain a small number of gluons, and "almost black" states (the Yin sector). 



Q I The excitations above the Yang vacuum have natural interpretation in terms of gluons. Analogously 

1^ ' the excitations above the Yin vacuum have natural interpretation as "holes" in the black disk - points 

at which an incoming gluon does not scatter with unit probability. We discuss in detail the spectrum 

^ ■ of the "parton model approximation" to the KLWMIJ evolution introduced in our previous paper, 

j^ . and prove that it is explicitly selfdual. This allows us to find explicitly the counterpart hole states 

in this approximation. We also present an argument to the effect that the end point of the evolution 

for any physical state cannot be a "grey disk" but must necessarily be the "black disk" Yin vacuum 

state. Finally, we suggest an approximation scheme for including the Pomeron loop contribution to 

the evolution which requires only the solution of the JIMWLK/KLWMIJ Hamiltonian. 
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1. Introduction 

In this paper we continue our investigation of the correspondence between the Reggeon Field Theory 
(RFT)[0, §, m and the modern incarnation of the ideas of gluon saturation [Q, ^, |^, as encoded 
in the JIMWLK-KLWMIJ functional evolution equation |, |, 0, |n|, |l|. In @ we have pointed 
out that the JIMWLK-KLWMIJ Hamiltonian can be directly interpreted as the Hamiltonian of RFT 
as arising from the first principles QCD calculation. The basic quantum field theoretical degree 
of freedom in this Hamiltonian are the unitary matrix R{x), which is the single gluon scattering 
matrix. In this respect this RFT differs from the original Gribov's template considered in 0, 0] and 



also recently in |14| which deals directly with scattering of color neutral physical objects. It is rather 



more conceptually related to Lipatov's approach ||15| in which the basic entity is the color nonsinglet 
reggeized gluon fl^ . 

We have shown that the reggeized gluon as well as other colored reggeon like states arise as the 
eigenstates of the RFT Hamiltonian in the partonic type approximation which preserves the number 
of s-channel gluons throughout the energy evolution. In the same approximation we have found 
the BFKL Pomeron as well as other Pomeron like and Odderon eigenstates. Thus the parton like 
approximation to the RFT reproduces the known perturbative QCD exchanges as well as contains 
other exchanges which only appear if multiple scattering of a single projectile gluon with the target 



is allowed. These additional exchanges grow faster with energy but their coupling is suppressed by 
powers oi as^. 

The spectrum of the full KLWMIJ (or equivalently JIMWLK) hamiltonian is not known apart 
from two eigenstates which both have zero eigenvalue. We will call these two states "white" or 
Yang, and "black" or Yin. The white state is simply the QCD vacuum state - the field strength 
in this state is zero. The scattering matrix of any projectile on this state is unity and it does 
not evolve with energy. The "black" state corresponds to the black disc limit. The probability 
density to find any gluon field configuration in this state is independent of configuration, and the 
scattering matrix of any projectile on it vanishes. These two states are related to each other by 
the Dense-Dilute Duality (DDD) transformation discussed in [^]. The DDD transformation was 



shown in [O, 18 1 to be an exact symmetry of the full RFT hamiltonian which should include the 



Pomeron loop effects H, |g, ^ |12|, |2|, g, |4|, ^. Although the JIMWLK-KLWMIJ hamiltonian 



is not self dual, the remnant of this symmetry survives on the zero eigenvalue subspace. The full 
selfdual RFT Hamitonian (which we will refer to as H ) is not known yet, although a number of 
selfdual extensions of the JIMWLK evolution has been proposed and derived in some approximations 
|2|, 1^, H] (see however @). 



In the present paper our goal is different. Rather than study specific approximation schemes or 
specific models, we aim to discuss some general properties of the spectrum of RFT, specifically those 
properties which follow from the exact selfduality which the full Hamiltonian of RFT must possess 



mm- 



The structure of this paper is the following. In Section 2 we recall the basic equations of the high 
energy evolution in the JIMWLK-KLWMIJ framework. We also give a derivation of the general form 
of the evolution in terms of n- gluon emission amplitudes. This section sets the stage and notations 
for the subsequent discussion. The relation of the high energy evolution kernel to the multigluon 
emission amplitudes to our knowledge has not been discussed in the literature so far. It is useful for 
the understanding of the general structure of the evolution but is not essential for the rest of the 
discussion in this paper. 

In Section 3 we show that the Hilbert space of the full (and as yet unknown) RFT Hamiltonian 
splits into two orthogonal subspaces which we dub the Yin and the Yang Hilbert spaces. The Yang 
states are those close to the white vacuum. They contain a small number of gluons and have 5-matrix 
close to unity. The dual image of this space is the Yin space. These are states close to the black 
vacuum. To characterize these states we find it natural to introduce the concept of a "hole". The 
hole at a transverse coordinate x is a state close to black disk, such that a gluon that scatters on 
it has unit probability to scatter unless it hits precisely at the point x, in which case it has a finite 
probability to propagate through the target without scattering. This concept naturally arises in the 
context of JIMWLK evolution, where the single gluon scattering matrix S acts as a hole creation 
operator. 

We illustrate the Yin- Yang structure within the partonic approximation to the KLWMIJ evo- 



lution introduced in [13|. We show that this approximation is itself selfdual, albeit under duality 



^The parton like approximation is designed for the situation when the projectile is small, while the target is far from 
being black. Therefore, just like the original BFKL approximation it is not applicable at asymptotic energies, and in 
fact it violates the unitarity of the scattering amplitude. 



transformation which is somewhat modified relative to that discussed in [17|. Using the results of 



1 13] we find explicit eigenfunctions in the Yin sector with up to two holes. 

In Section 4 we discuss how the eigenfunctions of the evolution Hamiltonian enter the calculation 
of physical scattering amplitudes. To this end we recall the restrictions on allowed states of RFT 
imposed by their interpretation from the point of view of scattering of physical gluons. We show 
that these restrictions are such that any physical state must have a nonvanishing projection on 
the Yin vacuum state. We also write down the expression for the physical s - matrix in terms 
of the eigenfunctions of H . Considering carefully the restriction on physical states we argue 
that asymptotically at large energies all physical states become black. Thus if an additional "grey" 
vacuum state exists it decouples from the evolution of physical states. 

Finally in Section 5 we conclude by discussing an approximation scheme to H^^'^ which includes 
Pomeron loops but requires only the knowledge of the solutions of }{JI^wlk ^ 

2. The high energy evolution and the gluon emission amphtudes 

The process we are interested in is the scattering of a highly energetic left moving projectile consisting 
of gluons on a hadronic target. We are working in the lightcone gauge natural to the projectile wave 
function, A~ = 0. In this gauge the high energy scattering matrix of a single gluon projectile at 
transverse position x on the target is given by the eikonal factor ^ 

S{x;0,x~) = V expli f dy- T" a^{x,y-)\ ; 5(x) = 5(x; 0, 1) . (2.1) 

where T^ = if"-^'^ is the generator of the SU{N) group in the adjoint representation and xi - the 
transverse coordinate. 

The field ax is the large A~^ component created by the target color charge density. It obeys the 
classical equation of motion and is determined by the color charge density of the target Pt{x) via 

|3|, n, 0, lii 



a9r{x,x-)T<^ = g^ ^{x - y) [s\y-^,x-) pUy,x-)T- S{y-^,x-)] . (2.2) 

For a composite projectile which has some distribution of gluons in its wave function the eikonal 
factor can be written in the form analogous to S{x), see |]l] 



S^[aT] = [dpfW^'Ipp] expji/" dy- f cfx p1,ix,y~) aUx,y-)\ . (2.3) 

Here the subscripts P and T refer to the projectile and target resectively. The quantity pp{xi) is 
the color charge density in the projectile wave function at a given transverse position, while W^[p\ 
can be thought of as the weight functional which determines the probability density to find a given 
configuration of color charges in the projectile. For a single gluon p^{xi) = T"-5^{xi — x^), and eq.( |2.3| ) 
reduces to eq.(2.1). 



In our convention the variable x is rescaled to run from to 1. 



The total S-matrix of the scattering process at a given rapidity Y is given by 

S{Y) = I Da^ W^^[aT{x,x-)] J^^.YoIM^,^')] ■ (2-4) 

In Eq. ( |2.4| ) we have restored the rapidity variable and have chosen the frame where the target has 
rapidity Yq while the projectile carries the rest of the total rapidity Y. Lorentz invariance requires 
S to be independent of Yq [|^, |31| . 

The high energy evolution is given by the following general expression^: 

6 



-^S = j Da9,Wl[aT{x.x-)] H^^^ 



ax, 



5aT 



'^y-yMt{x,x )]. 



(2.5) 



Here H stands for a generic kernel of high energy evolution, which can be viewed as acting either 
to the right or to the left, as it should be Hermitian: 
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W'^iar] . (2.6) 



A complete quantitative treatment of the evolution is yet unavailable. It has been shown however 
in [^ that in order for the total 5-matrix to be Lorentz invariant and symmetric between the 
projectile and the target, the evolution kernel H^^^ must be self dual. That is it has to be invariant 
under the Dense-Dilute Duality (DDD) transformation 

A A 
a"(x,x-) ^ i-— — , — — ^ -ip%x,x-). (2.7) 



6p°-{x,x ) 



5a"'{x,x ) 



The derivation of the kernel is available in two limits - the high density limit where the color 
charge density in the wave function is assumed to be 0{l/a) and the low density limit where the 
number of gluons in the wave function is assumed to be 0(1). Both these limits break the symmetry 
between the target and the projectile [pS|]. 

In the limit when the color charge density of the target is parametrically large {p"" = 0{l/as)) 
the kernel is given by the JIMWLK expression |jlO|, O] 



H'^''^"^^ = ^ d'zQnz)Qt{z) = ^ i <fzQ:{z)Q:{z) 



(2i 



where the hermitian "amplitudes" Qi{z) are defined as 
Q^{z) = j d\^^^:^[S-\z) - S'^\x)]j'j,{x)- Qliz) 



[X - z 



d'x 



{x - z)i 
{x — z)"^ 



[S^''{z)-S'"'{x)]Jl{x). 

(2.9) 

The generator of the right color rotations Jr and the generator of the left color rotations Jl are 
defined aslE^l 



JUx 



'R 



-''{^W^°Ssfe 



Jl{x) = -iT{T''S{x) 



5S^{x) 



ji{x) = [s{x)jR{x)r. 

(2.10) 

^We have adopted the conventions in which the spectrum of the Hamiltonian H^^"^ is positive, and not negative as 
is standard in the recent hterature, see for example [jl3|. This aUows us to think about the system described by H^^^ 
as of a standard stable quantum mechanical system. 



In the representation where the JIMWLK Hamiltonian acts on functionals of a{x , x), the rotation 
operators are 

J«(x) = ^ , ^ -; J£(x) = , , ^ . (2.11) 

^^ ' 5a{x,x- = iy ^^ ' (5q(x,x" = 0) ^ ' 

The JIMWLK Hamiltonian is usually written in a form different from eq.( |2.8D , see for example [p^ ]. 
It is however straightforward to check that eq.( |2.8D is indeed equivalent to the standard JIMWLK 
kernel by using the unitarity of matrix S. 

The KLWMIJ Hamiltonian (the dilute target limit) can be written in a very similar manner 

jjKLWMij ^ ^ j d^zQ1{z)Q^dz) = ^ j <fzQt{z)Q'i{z) (2.12) 

where the hermitian "amplitudes" Q^^z) are defined as 

j2„{^ ~ ^)i iDab/x r>ab/\] Tb/\. r\a/ \ ^2^\^ ~ ^)i i jybar ^\ r>ba/\TTb, 



QUz) = J d'x^^^-^[R^\z) - R^\x)]rj,{xy, Q-{z) = J d'x^—^[R''^{z)-R'^x)]JUx) 

(2.13) 
where R {x) is the color charge shift operator 

ii(x) = Pexp{/dx-^;^T«}. (2.14) 

The generator of the right color rotations Jr and the generator of the left color rotations Jl are 
defined as 

J^(x) = -tr|i?(x)T'^^^|, Jl{x) = -tv!^T-R{x)j^Y Jl{x) = [R{x)jR{x)f, 

(2.15) 



or in the representation analogous to eq.(2.11) 



J«(x) =/(x,a;- = 1); J£(x) = /(x,x" = 0) . (2.16) 

The JIMWLK and KLWMIJ Hamiltonians are clearly related by the DDD transformation. 

In the rest of this section we present a framework for derivation of these formulae which gives a 
more general perspective on the high energy evolution. 

Consider a wave function of a hadron at certain rapidity j^). We are interested in calculating 
expectation values of a variety of observables in this wave function. In particular we want to calculate 
average of the S'-matrix - an operator which depends only on the charge density operators p"'{x). In 
general we will consider any operator that depends on p: 0[p{x)]. As we have discussed in detail 



in [12, 22 1 such a calculation can be represented as a path integral over p{x,x ), where x is the 
ordering variable. 

(^|0[p(x)][^) = Dp{x,x-)W[p{x,x-]0[p{x,x-)] (2.17) 



with a weight functional W, where the function 0[p{x,x~)] is defined by expanding the operator O 
in power series in p and assigning the value of x~ to each factor of p in the order it appears in the 
product. The exact value of x~ does not matter as the correlators p{xi)...p{x~) depend only on the 
ordering of x^ but not on their exact values. 



It is convenient to represent this expectation value in a different way. Define the charge density 
shift operator 



R 



5(3 



Pexp<^ / dx" / (fxp'^{x) 



(2.18) 



6l3''{x',x) 

where the path ordering is along the x~ direction and the integration is over the entire x~ axis. We 
can then write 



d[p{x)] -- 


= R 


' 5' 


0[/3(x,x-)]|^= 


=0 ; 


where 




r 


" 5 ' 


= (^1^ 


" 6 ' 



(*|6[/5(x)]|^) 
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1^) 



Dp{x,x )W[p{x,x ]R 



O[Pix,x-)]\p=0 (2.19) 
^ ^ (2.20) 



6p 



Therefore as far as the averaging over the wave function is concerned one always has to calculate the 
average of the operator R rather than that of an arbitrary operator O. Curiously, this expression 
can be written in the form very similar to that for the S'-matrix derived in [p^ j. The wave function 
^ in the Fock basis is simply a collection of creation operators acting on the vacuum 



1^) = ^[aT] |0) . 



(2.21) 



The operator R when acting on the gluon creation operator at transverse coordinate x rotates it by 

the "c- number" matrix 

( r A 1 

(2.22) 



R^ix)=Vexpydx-T- — 



x ,x 



Thus 



r 



6(3 



{^[a]\^[Rpa^). 



(2.23) 



Consider now the same hadron boosted to a slightly higher energy so that the rapidity of all the 
partons changes by AY. The wave function of the hadron changes. This change can be represented as 
the action of the unitary "cloud operator" C[p,aA,ci/^] where by a a and a]^ we cumulatively denote 
the annihilation and creation operators of gluons at rapidities in the new part of phase space opened 
up by the boost. 

1^) -^ C|^). (2.24) 



The explicit form of the operator C is known in the low density regime (see |^4[ for a detailed 
discussion). Our derivation here is general, not restricted to the dilute limit and does not require 
this knowledge. The state |^) does not contain the soft gluons created by a^: 



OA 1^) = 0. 
The color charge density changes due to the charge of the extra emitted gluons by 

Ap^ix) = a^(x)r'^aA(x). 
The density shift operator becomes 



Ra 



6(3 



Pexp<^ / dx' d^x[p"-{x) + Ap"(x) 



(5/3"(x-,x) 



(2.25) 



(2.26) 



(2.27) 



The expectation value of the same observable in the boosted wave function is now given by 

6 



(*lC^O[/(x) + A/(x)]C|^) = Ta 
6 



6p 



Ta 



<5/3 



{^\C^Ra 



513 



O[/3(x,x-)]|;3=0 
C\^). 



Note that by virtue of eq.(2.25) 



Ra 



sp 



I^) = R 



513 



1^) = *[i?/3a^]|0) 



The operator T therefore changes as a result of the evolution as 

Ar = Ta - r = (0|^[a]| C^^ACi?^ - l W[i?a^]|0). 



(2.28) 
(2.29) 

(2.30) 
(2.31) 



£.t 



The operator i?A C R^ has all the charge density operators and field operators rotated by Rp relative 
to the operator C 

-RaC[/5(x), aA{x), a}^{x)]R^^ = C[Ri3{x)p, Rf3{x)aA{x), i?/3(x)a^(x)] . (2.32) 

Also, as shown in |£2[, every operator p inside the averaging, when acting to the left /right, becomes 
the operator of the left/right rotation JL,R{R(i) ( p.l6 ), plus the ordering of the operators gets reversed. 
Let us formally expand the operator C in the Fock basis of the operators a a 



C[p{x), aA{x), a'^{x)]\OA) = ^ c„[/5] a^^xi) ■■■ a^' '^(x„)|Oa) 



(2.33) 



where |0a) denotes the Fock vacuum of the operators a a- Here the index n denotes the number of 
soft gluons as well as their color indices and transverse coordinates, so that less concisely 



ai,a2...,an 



-t [Al 



J)l,b2...,bn 
■^Xl,X2,...,X„ 



Then we have 



AF = 1^ 4[JL]i?J^"^(xi) ••• i?J"""(rc„)c„[Jfi] - ll F 



(2.34) 



(2.35) 



where the color indices of R[^ are contracted with those of the coefficient functions c. In the dilute 



limit only one gluon is emitted per one step of the evolution |12] and the KLWMIJ equation is derived 
such that only the term n = 1 is kept in Eq. ( |2.35| ) . In general however it is clear that multiple gluons 
are emitted, even though the explicit form of C is not known in a general situation. 

We can obtain another representation of AF concentrating on Hermitian observables O only. For 
these observables 



2F 



5(3 



{^\R 
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+ R^ 



513 



m 



(2.36) 



For any unitary A one has 2 — A — A^ = {1 — A^){1 — A) and we can use this "optical theorem" to 
rewrite eq.(2.31|) 



AF = -^(0|^[a] (l - C^RaCrI^ {l - RaC^ R^c') ^[^/3«1 



l^QlQnT. (2.37) 



The second equality is obtained by inserting the resolution of identity on the soft gluon Fock space, 
1 = ^„ \n^) (tt-aI with the n soft gluon state defined as 

|nA) = a^ \xi) ■ ■ ■ a^" t(x„) |0a) (2.38) 

and 

Qn[Rp] = (nA|l - C^jR,aA,ai] C[Jl, i?^ aA,i?/3 aL] |0a) • (2.39) 

Each Qn depends on transverse coordinates and has a set of n Lorentz and color indices which we 
do not indicate explicitly. 

Since r[^] by definition is the Fourier transform of W[p\ with p and 6/6P being conjugates, 
eq.( |2.20| ) , the evolution of W[p] is obtained from the evolution of r[jg] simply by replacing /3 — *• p 
in the evolution kernel. The infinitesimal form of the evolution is obtained by expanding Qn to the 
lowest order in the evolution parameter Ay. Thus the "weight function" W[p] evolves according to 

-^W[p] = H^^'^Wip] (2.40) 

with 

H''''^ = ^ ^ E QnlR] Qn[R]\AY=0 (2.41) 

n 

where the unitary operators R{x) is defined in terms of the derivatives with respect to the color 
charge density eq. (|2.14 ). 



To construct H explicitly we need to know the form of the cloud operator C. However there 
are several properties of the evolution which are independent of explicit form of C. First, since the 
evolution operator is the sum of squares, it is clearly positive definite^. Second, this form of the 
evolution explicitly preserves the signature symmetry R — > R^ as defined in [O] irrespective of the 
explicit form of C, since hermitian conjugation of Qn simply takes R into its transposed. 

The meaning of Qn is the amplitude for creation of n gluons as a result of the evolution to higher 
rapidity. In the KLWMIJ limit ( |2.12D only one gluon is emitted in one step in rapidity, thus the single 
gluon emission amplitude determines the evolution as in eq.( ^.13 ). In this case the index n denotes 



the color, the transverse polarization and the transverse position of the emitted gluon. The form of 
the emission amplitude eq.( p.l3 ) is precisely the same as that of the inclusive single gluon emission 



amplitude for the scattering on the target described by the scattering matrix R |35]. Beyond the 
KLWMIJ limit more than one gluon is emitted in one step of the evolution. For example in order 
to obtain the JIMWLK limit an arbitrary number of gluons would have to be emitted in one step of 



the evolution. This is consistent with the cloud operator being a Bogoliubov type operator |36]. 
3. The Yin and the Yang 



Our main point of interest is how to understand the evolution of the physical S'-matrix eq.( |2.4D in 
terms of solutions of H . The natural strategy for calculating the S-matrix is the following. First 



4n 



The derivative in the definition eq.(2.41) does not affect the conclusion about positive definiteness, since the product 



QnQn vanishes at AY — while it is positive for any nonzero AY. 



find eigenfunctions of the RFT Hamiltonian, 

^iJFT|^«,^ = uj,\¥i) (3.1) 

wliere tlie index li labels degenerate eigenstates corresponding to the eigenvalue Wj. The evolution 
of each eigenfunction is given by 

|$(»)y = e^'^'^l^'*). (3.2) 

Since H^^^ must be hermitian, its eigenfunctions must form a complete orthonormal basis. Next 
expand S at initial rapidity in this basis 

and similarly expand W^: 

^^ = E 7.7* (*ll (3-4) 

The total cross section at rapidity Y is 

where e-^ ' is the overlap matrix between degenerate eigenstates at a given value of Wj. 

In this discussion as well as in the rest of this section the scalar product and the norm of the 
eigenfunctions |^j) are defined in the standard way 

($,|M/,,) = j Dp^l[p]^,[p]- (3.6) 

Note that S and W are not necessarily normalized to unity with respect to this norm, thus ^^ 7/^ * ^f 7^ 



1 etc. The general properties of W and S have been discussed to some extent in |13]. We will return 
to this question in the next section. 

The question we are addressing in this section is how the DDD symmetry affects the spectrum 
of RFT. Since RFT is formulated as a regular quantum field theory with a hermitian Hamiltonian, 
quite generally any symmetry of this Hamiltonian must be represented in the spectrum in either 
Wigner or Goldstone mode. That is it is either the symmetry of the spectrum, or is spontaneously 
broken in the vacuum. In the latter case the vacuum has to be degenerate and the Hilbert space has 
to split into direct sum of two Hilbert spaces each one built over the appropriate vacuum state. 

Is DDD spontaneously broken or not in RFT? A little thought is required to see that the answer 
is affirmative. The simple argument is the following. The QCD vacuum state does not contain gluons 
and is Lorentz invariant - therefore it does not contain gluons in any frame. This state therefore is 
a vacuum of the RFT - the S-matrix of scattering it on any hadron is always unity and does not 
depend on rapidity {ojq = 0). This is the white "Yang" vacuum of RFT. On the other hand the DDD 
transformation as shown in [O] does not leave this state invariant, but rather transforms it into the 
completely black state - the "Yin" vacuum. Common sense of course tells us that the totally black 
state is also invariant under the high energy evolution, since the probability of any hadron to scatter 
on it is strictly equal to unity. Beyond the simple common sense the degeneracy of the Yin and Yang 
vacua is ensured by the fact that DDD is the symmetry of the RFT Hamiltonian. 



3.1 The two vacua 

The two vacuum structure is present already in the KLWMIJ and JIMWLK Hamiltonians, even 
though the DDD transformation is not a symmetry of each one of them separately. Consider for 
example the KLWMIJ Hamiltonian. It has been derived in the limit of a dilute system and thus 
must by definition contain the Yang vacuum. This is indeed the case. The eigenstates satisfy the 
Schroedinger type equation 

jjKLWMlJ^p-^ = a;^[p]. (3.7) 

The form of ff^^WJ^'^iJ ^s sum of squares ensures its positive definiteness. Thus the only zero 
eigenstates are those annihilated by all Qi in eq.( 2.13|) . We find useful in the following to use both 



the "configuration space basis" and the "momentum space basis". The first one is the basis of 
eigenstates of (mutually commuting) operators p'^{x,x~): 

p^ix, X-) Iq^ix, x')) = q^ix, x~) |g"(x, x')) (3.8) 

and the second is the basis of the eigenstates of the functional derivative operators 

— — |p'^(x, x~)) = p'^ix, x~) \p^{x, x-)) (3.9) 

p"-[x,x ) 

which obviously satisfy 

(g"|p") = expli dx^d'^xq''{x,x~)p''{x,x-)\ . (3.10) 

The Yang vacuum does not contain gluons and is thus annihilated by the color charge density 
operators. This obviously corresponds to the state^ \Yang) such that 

{q\Yang) = 5[q{x,x^)\; {p\Yang) = 1. (3-11) 

This state is indeed annihilated by the action of J%{x) and J£(x), and therefore also by the action 
of Qf (x) 

J^(^)(x) \Yang) = Q^{x) \Yang) = 0. (3.12) 

This is however not the only state annihilated by all Q^{x). The other such state is an eigenstate of 
all R{x) with the eigenvalue independent of x. We can choose this eigenvalue to be equal to one 

R''^{x)\Yin) = d^^lYin). (3.13) 

Since the operator R depends only on the derivatives with respect to p, the state \Yin) has the wave 
function which is independent of p 

{q\Yin) = 1; {p\Yin) = 6[p{x,x~)] . (3.14) 

The physics of \Yin) is precisely that of the black disk. The wavefunction does not depend on p, and 
thus all configurations of charge density are equally probable. Note that the actual derivation of the 



^For the purpose of the discussion in this section the overall normalization of the wavefunctions is not important. 
We will come back to this question in the next section. 
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KLWMIJ limit assumes that the charge density is smaU, thus strictly speaking \Yin) is outside the 
range of validity of the approximation. Nevertheless this state is a vacuum state of jj^i^^^^^J which 
appears as the remnant of the DDD symmetry of H^^^ . 

Considering the JIMWLK hamiltonian we find a very similar situation. }{JJMWLK j^g^g ^.^^ 
vacua, which are in fact the same as the vacua of jj^^^J^^J _ One state is annihilated by the 
dual rotation operators J'fi(i^)ix). The wave function of this state does not depend on a{x,x~) and 
therefore by eq.( |2.2| ) also on p"(x,x^). This state is nothing but the black disk \Yin): 

J^(^)(x) \Yin) = Q^{x) \Yin) = 0. (3.15) 



The existence of this state in the JIMWLK framework was first pointed out by Weigert in ||10|] . The 
other vacuum state is the eigenstate of S{x) with eigenvalue one. Since S is the scattering matrix 
of a single gluon on the wavefunction in question, this state clearly corresponds to physical vacuum 
and up to a possible normalization constant is equal to \Yang). 



S'^'ix) \Yang) = 5"" \Yang) . (3.16) 

The JIMWLK Hamiltonian was derived under the assumption of large color charge density, therefore 
strictly speaking \Yang) is outside its range of validity. Nevertheless in the framework of JIMWLK 
evolution this state as well as states close to it (with small charge density) have been studied in the 
original papers [|l^] as well as more recently in relation to the Odderon exchange |37]. 



3.2 The doubling in the spectrum 

As noted above, the KLWMIJ and JIMWLK Hamiltonians separately are not self dual, thus one does 
not expect exact degeneracy in the spectrum of either. Nevertheless, the spectrum of each contains 
two towers of states built above the two vacua. To illustrate this point we will discuss the KLWMIJ 
evolution. The exact spectrum of the KLWMIJ Hamiltonian is not known. However in |13] we have 
studied it in the partonic approximation. This approximation is obtained by expanding R around 
unity while preserving the ultraviolet and infrared finiteness of the KLWMIJ Hamiltonian in the 
leading order approximation. The leading order Hamiltonian is given by 

with 

Qr'^iz) = Jd'x^^tr ([T^i?(x) - R{z)] ^^) (3.18) 

with R = R — 1. In this approximation we have found a "vacuum state" of the form eq.( [4.lD with 
E[i2] = 1. This is just \Yang). We have also found eigenstates of the form 

G^= I e^'^^r,^,,R^\x)\Yang) (3.19) 

J X 

where Vj^^ is a color projector onto a representation labeled by A, and 

G\= P^^^ I R''\u)R^'^{v)¥g{u,v)\Yang) (3.20) 
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with P^^jj - a projector from the product of two adjoint representations onto a representation labeled 
by i. The wavefunction ^Uu,v) for any i satisfies the BFKL equation. 

As explained in detail in |13], from the t-channel point of view each factor of t] R corresponds 



to a t - channel exchange with the color quantum number in representation A. Each wavefunction 
in eq.( 3.19 ) describes one such exchange, whereas the wavefunctions in eq.( p.20 ) describe a pair of 



such exchanges which include the BFKL pomeron and odderon trajectories. 

On the other hand when interpreted in terms of the s - channel gluons, every power of R in the 
wave function of RFT corresponds to a gluon in the QCD state. For example the state 

|G"^(x)) - R''^{x)\Yang) (3.21) 

corresponds to the QCD state with a single gluon at the transverse position x which changes its color 
index from a to 6 in the process of scattering. Analogously 

|G"'''^'^(x,y)) ~ R^\x)R''''{y)\Yang) (3.22) 

is the RFT representation of the QCD state with two gluons at transverse positions x and y. Thus 
the wavefunctions of eqs.( p.l9| , 3.20| ) contribute to the evolution of the scattering amplitude of a one 



and two gluon states respectively. These states therefore naturally belong to the sector of the Hilbert 
space which is close to the Yang vacuum. 

Since ifP"''* is homogeneous in R one may expect to get from it some information about the Yin 
part of the spectrum also. This is indeed the case. First of all it is obvious that the Yin vacuum is 
also the vacuum of i/P°-^*. More importantly it turns out that ifP"^* is invariant under a modified 
version of the DDD. In particular it is a matter of simple algebra to check that //p"''* is invariant 
under the following transformation 

R{x) ^ [d'y^{x,y)-I—- .^ ^ d'R(^x) (3.23) 

with -m{x, y) = 2tt ^^[{^ ~ y)^/^'^]- This is the linearized version of the DDD transformation. 
Thus we immediately infer that the states in the Yin sector of the form 



and 



Gl = m [ ln[(n - x) V'] ln[(^ - y) V'] ^^(x, y) -^^ -J— \Yzn) (3.25) 



SR^^iu) 5R^<^[v) 



are eigenstates of ifP"^*. The corresponding eigenvalues are identical to those of the corresponding 
states in the Yang sector. 

We note that although the parton like approximation is not identical to the BFKL approximation 
to the KLWMIJ evolution, the invariance under the modified duality transformation is shared by the 
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BFKL Hamiltonian. The BFKL approximation is equivalent to expanding the factors R{x) in the 
KLWMIJ Hamiltonian to lowest order in 6/6p. The resulting Hamiltonian is^ 

^BFKL ^ ^ /'^2^gBFKLa(^)gBFKLa(^) (326) 



with 



gBFKLa(^) = I d^:, i^^ f-bc 



(x-z) 



2 



Sp'^ix) 5p^{z) 



p%x) (3.27) 



This expression is invariant under the duality transformation 



^ '' cPyl-{x,y)p%y)- p«(^)^a2— ^. (3.28) 



Thus the exact doubling of solutions corresponding to the Yin and Yang sectors exists also in the 
BFKL approximation. 

3.3 Holes in the black disk: the JIMWLK perspective 

The structure of two vacua and two "towers" of excitations is also supported by the JIMWLK 
limit. This is inevitable since the JIMWLK and KLWMIJ spectra are unitarily equivalent as the 
two Hamiltonians are related by the unitary DDD transformation. It is nevertheless instructive to 
analyze the spectrum from the JIMWLK viewpoint. 

First, as discussed above the two vacua of JIMWLK are exactly same states as the two vacua 
of KLWMIJ. This time however, it is the black disk state - the \Yin) which is the natural vacuum 
of the JIMWLK evolution, while the presence of \Yang) is "accidental" and is the remnant of the 
DDD symmetry of H^^^ . 

To get an idea about the nature of the eigenstates we can employ the approximation dual to the 
parton model eq.( |3.17|) . 

with 

Qr\z) = Jd'x^^^tr (^[T^,S{x) - S{z)] ^) (3.30) 

with S = S — 1. The eigenstates of this Hamiltonian are obviously DDD duals of the "reggeon" and 
"pomeron" states found in [13|. The states of the Yin "tower" are 

H^ = f e''^''r]^aS^'^ix)\Yin) (3.31) 

J X 

Hl= P^:i f S'^\u)S^\v)%{u,v)\Y^n) 

Ju,v 

^Though the BFKL Hamiltonian can be formally written as the sum of squares, as is well known, its lowest eigenvalue 
(in our convention) is negative. The resolution of the paradox is in the fact that the BFKL eigenfunctions are not 
normalizable and therefore positive definiteness of the Hamiltonian does not hold. 
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and the like, while the states in the Yang tower are 

H^ = I e''" Vcd j^ \yang) (3.32) 



y 



Hl= n^ I ln[(u-x)V]ln[(«-y)V]n(^'y)T^i7TTa^l^«^5) (3.33) 



'x,y,u,v 



etc. 



The interpretation of the states of eg. ( p. 31 ) is quite amusing. Consider for example the dual of 
a single gluon state eq.( |3.21| ). 

\H''''{x)) = 5'^*(x) \Yin) . (3.34) 

This state is black almost everywhere except at the transverse position x. To see this imagine that 
we scatter on this state {W = \H"''^{x))) a single gluon projectile at transverse position z. The 
projectile averaged single gluon S-matrix is simply T,[S] = S'^'^{z). To calculate the scattering matrix 
for such a process we have to integrate this S with the weight function W corresponding to the state 
eq.(pl) 7 

S = j DSS'^\z)S'^\x) = ^_^5,_,5-5^<^ (3.35) 

with the "Kronecker delta": 5q = \. Thus multiplication of the black disk state \Yin) by a factor 
S{x) produces the state which is not completely black for a gluon that impinges on it at the point x. 
It is natural to think of such a state as a hole in the black, see Fig. 1. We will adopt this terminology 




Figure 1: The Yin-Yang structure: duality between gluon excitations above the white vacuum and the hole 
excitations above the black vacuum. 

and will call S{x) - a hole creation operator. Thus the excitations in the \Yin) sector can be though 
of as the "reggeized holes" and their bound states - the "hole Pomerons". Even though formally 
from the point of view of the eigenstates of H there is complete symmetry between the gluon 
and the hole states, physically of course they are very different and they contribute quite differently 
to physical amplitudes. We now turn to the discussion of this question. 

4. The RFT Hilbert space and the physical states 



To calculate physical amplitudes using eq.(3.5) we need to understand how to expand a weight 



functional W that describes a physical s tate in the basis of eigenfunctions of H . As mentioned 



^Eq.(3.34) in fact defines an SU{N) ® SU{N) multiplet of states corresponding to the indices ab. The same is true 
for the projectile averaged S'-matrix . For scattering of a gluon with a given color on a given state we should specify all 
four indices appropriately. We will however keep the indices arbitrary as it does not cause any additional complications. 
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above the functional W cannot be arbitrary but rather must satisfy some conditions. This issue has 



been discussed before in |13] and |2^. We repeat the main points here since they are important for 
our subsequent discussion. 



4.1 Normalization of physical states 



First, as shown in |29] in order for the correlators of powers of p'^ to be consistent with the quantum 
commutation relations of the operators p" the functionals W must be of the form 



W[p] = ^R]6[p{x,x')] 



(4.1) 



with an arbitrary functional S which depends only on the charge shift operator or the "dual Wilson 
line" R. Eq.(4.1) is a restriction on the Fourier transform of W. Inverting it we find 



/«,exp{i/..-AA„-).»,.,:,n}H'H^E|5, 



(4.2) 



with S defined as in eq. (|2.l| ). Thus the Fourier transform of W must depend only on the Wilson 
line S and not on any other combination of the conjugate variables a. This has a transparent 
physical interpretation. The functional Fourier transform of W is precisely the projectile averaged 



scattering matrix for scattering on the "external field" a. The restriction eq.(4J) simply means that 
this scattering matrix must be a function of scattering matrices of individual gluons constituting the 
projectile, and does not depend on any other property of the external field. The normalization of 
the functional W is determined by requiring that for a = the scattering matrix is equal to unity. 
This leads to the overall normalization condition 



S[i2 = 1] = 1, 



j DpW[p] =1 



(4.3) 



Further general properties of S follow from its identification as the projectile averaged scattering 
matrix. Suppose for example the wave function of the projectile at the transverse position x has the 
form (we take there to be exactly one gluon at the point x) 



N^-l 



\^x) = Y^ Ca{x)\a,x) 



(4.4) 



a=l 



The scattering matrix operator when acting on ^ multiplies the gluon wave function by the matrix 



S 1 26]. We thus have 



Af2-1 



(^,|5|^,) = ^ Ca{x)Cl{x) S^\x) 



(4.5) 



a,b=l 



This generalizes to states with more than one gluon. Expanding the functional S in Taylor series in 
R we obtain the general form compatible with its interpretation as the projectile averaged scattering 
matrix 



oo oo 



^i-^J — 2.^ Z-^ K'{ai}...{a'"}(^l) •••)^m)C^{fti}...{bm}(xi, ..., j;m) 



m=0 ni=l 



n: 



i=l 



R^l''UxA...R 



at . bl 



(x.) (4.6) 
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where {a^} denotes the set a\,al, ...,a^^, etc. Simply stated this is the 5-matrix of a state whose 
quantum mechanical amplitude to have rij gluons with color indices a^^, ..., a^, at the transverse coor- 
dinate Xi with i = 1, ..., m is Cr^i j r^m|(xi, ..., Xm)- Thus for example the coefficients of "diagonal" 
terms with a*- = 6*- must be positive, since they have the meaning of probabilities to find the con- 
figuration with particular color indices in the projectile wave function. One should note that in the 
derivation of eq.( |4.6D we have assumed that the projectile is in a pure state and is described by a wave 
function. If instead the projectile is described by a mixed density matrix a wider set of functions S 
is allowed. We will not enter into a more detailed discussion here as we will not have to use these 
conditions below. 

As is obvious from the discussion above, the normalization conditions on physical states are very 
different than the conventional normalization of states in the Hilbert space of RFT eq. (|3.6| ). 

Next consider the eigenvalue equation of the JIMWLK Hamiltonian 

hJIMWLK^^^ ^^^^_ (4.7) 

Since //-'^^M/Lk ^g proportional to the generator of rotations, integrating this over the group with 
any group invariant measure leads to 



/ 



Dp^i[S] = 0, for w» y^ (4.8) 

Another property of eq.( [4.7D is that its left hand side vanishes at S = 1. Thus we also have 

^,[S = 1] = 0, for io^ ^ 0. (4.9) 

Since the eigenfunctions of }{^^^^^^-J are functional Fourier transforms of the eigenfunctions of 
jjjiMWLK ^ and since eqs.( [4.8D and ( [4.9]) are Fourier transforms of each other, we conclude that 
the eigenfunctions of ff^LWMU gj^^re these properties. This can of course be seen by considering 
directly the eigenvalue equation for the KLWMIJ eigenfunctions. 

Eqs.(4.8,4^ also hold for the wavefunctions of the complete H . Eq.( [4.8D is simply the 



statement that all eigenstates with nonvanishing eigenvalues are orthogonal to the Yin vacuum, 
while eq.([4.9|) is the orthogonality condition of these states with the Yang vacuum. The property of 
orthogonality of eigenstates corresponding to different eigenvalues has to hold for eigenfunctions of 
any hermitian operator, and thus also H^^^ . 

Eq.( [4.8D in conjunction with the normalization of the physical wave functions eq.( ^.3D has an 
immediate consequence that any physical state must have a nonvanishing projection onto a state 
with zero energy. We know that H has two zero energy states - Yin and Yang. We will assume 
for the purpose of the current discussion that those are the only two states with zero energy. We will 
discuss the possibility of existence of additional zero energy states in the next subsection. 

We now have to discuss in more detail the normalization of the Yin and Yang eigenfunctions. To 
write the properly normalized wave functions in a sensible way we have to introduce the infrared and 
ultraviolet cutoff on the system - the total volume V and the closest distance 1/A in the transverse 
plane. We also have to introduce the "infrared" and "ultraviolet" cutoffs in the field space. In 
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particular we will cutoff the p integrals at some large value v and at a small value A. With these 



cutoffs in place we can revise eq.( |3.11| , p.l4D appropriately 



Go = {q\Yang) = H, A^/^ 5(g(x)); 
i/o = {q\Yin) =Ii^v-^'^. 



(4.10) 
(4.11) 



Since we would like to keep the DDD explicit at finite cutoff, the two cutoffs A and v are not 
independent. In particular, we want the state Hq to be the Fourier transform of the state Gq. This 
requires A = v~^ . Thus for the normalized eigenstates we have 



j DpGl = j DpHl = 1; J DpGo = v-""/^; fDpHo = v 

where n = VA is the total number of points in the transverse plane. 



(4.12) 



With this normalization the expansion of a physical state in the eigenfunctions of H has the 
form^ 

W[S] = 70^2 Go + mv-'^'^H^ + W2 Y. 7,G,[5] + t;-"/^ J] r^,Hi[S] (4.13) 



i:uJi>Q 



i: uJi>0 



with 



70 + r/o = 1 



(4.14) 



which follows from the overall normalization ( [4.3[) . The states Gj and Hi are "gluon" and "hole" 
exited states belonging to the Yang and Yin Hilbert spaces respectively. Assuming that the nor- 
malization of other states in the Yang tower Gj is similar to that of Go as far as the volume and 
ultraviolet cutoff is concerned and the same for Hi relative to Hq, we expect all coefficients 7^ and rji 
to be numbers of order unity. For S, which is the Fourier transform of W the expansion is analogously 

S[5] =7oW2i/o + r?ot;-"/2Go + W2 ^ ^i Hi[S] + y-''/^ Yl ^iGi[S]. (4.15) 

i:LJi>0 i:u>i>0 

Note that although the Yin and Yang sectors form separate Hilbert spaces in the thermodynamic 
limit, at finite value of the regulators the overlap between the states in the two sectors is finite. In 
particular 

{Yin\Yang) = f DpGoHo = u"" . (4.16) 

Thus for scattering of a projectile P on a target T we have 
S{Y) = (7o^ + v^) ill + ri) 



v^ - v-^ P r 

- % ^0 + 



V'' 



(4,17) 






u),Y 



ifr^ + riflf + [v-lflf + v-^njr^f) j DpG.H, 



^Compared to Eq. (3.4) we explicitly take into account the doubling of the spectra. Namely the index U runs up to 
two: eluons and holes. 
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or in the infinite volume limit 

i:uji>0 



ifvl + Vfll + 7f 7>" / DpG.H, 



(4.18) 



This formula tells us that if neither the projectile nor the target have an overlap with the white 
vacuum {-jj = -y^ = 0; rj^ = rj^ = 1) , the S - matrix at large rapidity vanishes, as only the positive 
energy eigenstates contribute. If either the projectile or the target have a white component, this 
means that ?/o t/q ^ 1, and therefore the S matrix tends to a constant value as y ^ oo. This is 
the situation when the projectile (or the target) has a vacuum component in its wave function - 
this component remains transparent even at asymptotically high rapidity. This of course is not the 
situation one usually is interested in. It is in fact easy to see that physical states that contain finite 
number of gluons are all orthogonal to the \Yang) state. For example a single gluon state (eq.( |3.2l"l )) 
properly normalized to satisfy eg. ([4. 3D 



{YanglC^x)) = v""^^ {Yang\R''\x)\Yang) = (4.19) 

by virtue of \Yang) being invariant under the SU{N) x SU{N) group that rotates the left and right 
indices of R. The same argument applies to any state which contains a nonzero number of gluons at 
some transverse position. On the other hand all these states have the same overlap with the \Yin) 
vacuum 

{Yin\G''''{xi)...G'"^{xn)) = v""^^ {Yin\R''''{xi)...R'"^{xn)\Yang) = (4.20) 

= v''/^5''''...6^'^{Yin\Yang) = y-''/'^ 6''\.r'^ . 

Thus for all physical states of interest 7o = and r/o = 1- 

4.2 Tertium non data 

The existence of the white and black vacua of H is not in doubt. An interesting question is 
whether the theory could have an additional "vacuum" a; = state. Such a state would be selfdual 
and presumably "grey" in the sense that it would lead to S matrix not equal to zero or unity. 

At first sight the existence of such a state is suggested by the gluon-hole duality itself. Indeed, 
we know that the state which contains initially a small number of gluons becomes darker as the result 
of the evolution. This darkening is illustrated for example by the growth of the number of gluons. 
The gluon number operator in the framework of the BFKL evolution is simply related to the Fourier 
transform of the charge density operator 

Ng{k)^^{p%k)p%-k)). (4.21) 

Within KLWMIJ evolution the charge density correlator satisfies the BFKL equation and thus grows 
exponentially with rapidity. The number of gluons therefore grows exponentially fast and the state 
becomes blacker. 
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What happens if instead one starts with an initial state close to the black disk limit and follows 
the evolution of the number of holes? The observable corresponding to the number of holes clearly 
is just the dual of the number of gluons 

The rapidity evolution close to the black disk limit is given by the JIMWLK equation. This leads 



to a linear BFKL evolution of the correlator in eg. ( 4.22 ) in a manner totally analogous to the linear 
evolution of the p correlators in KLWMIJ. So we conclude that close to black disk the number of 
holes exponentially grows with the BFKL exponent. 

It is tempting to conclude from this argument that the black disk is an unstable state and if one 
chooses an initial condition sufficiently close to it, this initial state will flow away from the black disk 
in the course of the evolution. The final destination for such a state then may be some state other 
than either black or white vacuum - a "grey vacuum". If such a state exists, one would expect its 
properties to be quite interesting. First, it would be invariant under the DDD transformation. The 
gluon operator is equal to unity when acting on the black disk state, eq. (|3.13| ). In this sense the black 
disk can be thought of as a condensate of gluons, which breaks the SU{N) ® SU{N) gluon rotation 
symmetry down to the diagonal SUv{N). The white state on the other hand according to eq.( p.l6 ) is 



the condensate of holes and breaks the hole rotation symmetry in the same fashion. One could expect 
the grey state to be selfdual, break no symmetries and have neither gluon nor hole condensates. In 
physical terms that would mean that the scattering probability of say a dipole of size r will be equal 
neither to zero nor to unity for any finite size r. One would of course still expect the probability to 
approach unity for r — > oo and zero for r ^ 0, but only asymptotically. The most natural realization 
of such properties would be if the state was conformally invariant and yielded correlation functions 
which behaved as a power of the distance. We note that these properties are consistent with the 
asymptotic behavior suggested in |^5|. Such a selfdual state also exists on the classical level in 
Gribov's Reggeon field theory [||, |2|, |lj], although its fate beyond classical approximation is unclear^. 

The hole proliferation argument however cannot be interpreted directly as pointing to the exis- 
tence of an additional vacuum. The hole states crucially differ from gluons in that they do not satisfy 
the conditions of physicality. In particular any state with finite number of holes has a zero norm in 
the sense of eq.( [4.3D , since it is by definition orthogonal to the \Yin) vacuum. Since the condition of 
eg. ( [4. 3D expresses the normalization of the QCD wave function of the projectile (or target), zero norm 
states necessarily contain negative probabilities and thus are clearly unphysical. A complementary 
feature of the hole states is that they are not orthogonal to the \Yang) state. Thus even if physical 
interpretation were possible, the QCD state involved would contain a vacuum component in its wave 
function. It would not be therefore of the type of state which one expects to get black as the result 
of high energy evolution. The proliferation of holes may be simply a manifestation of the fact that 
the "white" component in the hole wavefunction dominates the S matrix at high energies. 

Indeed, let us present an argument to the effect that the grey state if it exists does not play any 
role in the evolution of physical states. We have shown in the previous subsection that any physical 



®The conclusion of |l4| e.g. is that the asymptotic behavior is indeed grey but not related to the classical selfdual 
solution p — q = 0. 
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state, namely any state consisting of a finite number of gluons has unit projection on the \Yin) state, 
r/o = 1. It therefore follows from eq.( 4.14| ) that its projection onto any other state with unit physical 



norm defined by eq.([4.3|) vanishes^''. As regards the grey vacuum this means that either a) it does 
not exist, or b) it exists but its overlap with any physical state vanishes, or c) its physical norm 
vanishes, in which case it cannot be a final point for evolution of any physical state. Our present 
argument does not exclude the possibility c), namely that the finite point of the evolution is an 
admixture of the black vacuum and a zero physical norm " grey vacuum" . We view this situation as 
highly improbable, since in this case the "grey vacuum" itself is not a physical state, unlike the Yin 
and Yang states. We believe this possibility can be excluded by consideration of other restrictions 
on the nature of physical states as encoded in eq.( |4.6| ). With this proviso we conclude that a grey 
state even if it exists in the Hilbert space of H does not couple to evolution of physical states^^ . 



5. Discussion 

In this paper we have discussed the structure of the Hilbert space of H^^^ , the realization of the 
DDD symmetry and some general consequences of this for the high energy evolution of physical 
states. We have shown that the DDD symmetry is spontaneously broken leading to appearance 
of two zero "energy" states - the white "Yang" vacuum and the black disk "Yin" vacuum. The 
excitations in the Yang Hilbert space have the character of gluons (or rather reggeized gluons) and 
multigluon states. The dual of these states in the Yin Hilbert space can be thought of as holes in 
the black disk (or reggeized holes to complete the analogy). The gluon "creation operator" in the 
framework of the RFT is the color charge density shift operator R. It has the double meaning of 
creating the charge density of a gluon in the hadronic wave function and also as the scattering matrix 
of the gluon belonging to the projectile when it scatters on an external hadronic target. The hole 
operator S is the scattering matrix of a gluon impinging on the hadronic wave function. 

The gluon and hole operators are dual with respect to DDD and fit the standard pattern of 
order-disorder variables taken to the extreme. The operator R{x) creates an extra adjoint charge at 
the point x, while the operator S{y) measures the color field created by this charge. Perturbatively 
this field is a{y) (x gln{y — x)'^fj? where fj, is the infrared cutoff. Thus the phase of operator S{y) 
created by the action of R{x) is of order asln{y — x)'^fi'^. At distances |y — x| < Tc = l/zie^'"" 
this phase is very large. The operator R therefore disorders any configuration of S on the distance 



^''Eq.(4.14) assumes existence of two vacua only, but its generalization to a greater number of such states is straight- 
forward. 

^^This argument seemingly contradicts conclusions of [M] as well as those of the recent paper |3q|, which finds a grey 
state as the final state of the high energy evolution in a toy model without transverse dimensions. We note however 
that the properties of the evolution in the toy model of |B8] diff'er in a crucial way from those in QCD. In particular the 
Hamiltonian of [P8| is not hermitian, thus many of our arguments do not apply to this case. Also, and perhaps more 
importantly, as explained in ISSl the toy model Hamiltonian describes real processes of both emission and annihilation 
of dipoles in the wave function. In such a situation one indeed expects that the evolution should terminate on a grey 
state, since the dipole density in the wave function must literally saturate. In QCD on the other hand only emission 
of gluons is present but annihilation is not allowed, see pst] for discussion. Thus one does not expect saturation of the 
gluon density in the QCD wave function, although the rate of its growth is expected to slow down dramatically at high 
energy. 
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scale re- Since the scale re itself is proportional to the infrared cutoff, this is an extremely strong 
disordering effect. 

The Yin and Yang vacua fit naturally in the picture of order-disorder duality. The hole operator 
S is equal to unity when acting on the Yang vacuum, thus holes are condensed in \Yang). Likewise, 
\Yin) is the condensate of gluons. The symmetry breaking pattern in the two states is also different. 
The gluon operator transforms as an adjoint representation under the SUg{N) ^ SUg{N) transfor- 
mation as R —f Ug RUg. Although we do not know the full H Hamiltonian. we do know that 
this transformation is a symmetry of its KLWMIJ limit. Analogously the JIMWLK Hamiltonian 
is invariant under the SUh {N) ^ SUh {N) transformation S — > UhStlh. The vector subgroup of 
these two groups is the same, since under it both p^ and bjbp"" transform as adjoint representation. 
However other group elements are different. In particular in the Yin vacuum the gluon symmetry 
SUg{N)®SUg{N) is broken down to the diagonal SUv{N) but the hole symmetry SUh{N)(8)SUhiN) 
remains unbroken. In the Yang vacuum the situation is reversed. This situation is generic for theories 
involving order-disorder variables. Some symmetries which are represented a la Noether for one set 
of variables look topological when written in terms of the dual set, and vice versa. Although we have 
not found an explicit realization of say SUh{N) SUh{N) in terms of gluons, we may expect that 
the elements which do not belong to the vector subgroup have topological interpretation. 

The spontaneous breaking of these global symmetries should lead to the appearance of the 
Goldstone boson modes. Indeed in the framework of the parton like approximation to the KLWMIJ 
equation these Goldstone modes were identified with the reggeized gluons. By duality the reggeized 
holes are the Goldstone modes in the JIMWLK case. Whether this situation persists in the full RFT 
is at this point unclear. The full H is selfdual. Thus either both SU{N) (81 SU{N) groups are 
exact symmetries of the full RFT and the spontaneous breaking pattern in the vacua persists, or only 
the common diagonal subgroup remains the symmetry of H . In the latter case the Goldstone 
bosons only appear in the dense and dilute limits and are artifacts of these limiting cases. 

Although our discussion has been rather formal, we have been able to draw some interesting 
conclusions. In particular we were able to argue that the high energy asymptotics in the full H 



must be black and not grey. The formula eg. ( 4. 18 ) is also informative. For physical states (r/ = 1) 
the S matrix is given by 



SiY) = Y: e-^^ ^jfvT + Vflf + itlTv- JDpG^H, 



i: uji>0 



(5.1) 



One can infer from it for example that the approach of the S matrix to the black disk limit must be 
exponential. At large Y the leading contribution comes from the lowest eigenvalue that couples to 
the particular physical state. In principle in the case of a gapless spectrum this eigenvalue could be 
zero leading to a power like approach to the black disk limit. However if the only massless modes 
are the Goldstone bosons of the type discussed above, and therefore carry color quantum number, 
they will decouple from the S - matrix of any color neutral physical state. Their presence could 
nevertheless be felt in that the spectrum even at finite lo can be continuous. From this perspective 



we view the recent result about the asymptotics of the solutions of the Kovchegov equation [39| as 



very natural. It was shown in |39| that the approach of the S matrix of a dipole to the black disk 
limit does not follow a gaussian S oc e~^^ as suggested in [^], but is rather exponential with a 
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power prefactor, S oc Ye "^^ . From the perspective of eq. (|5.1| ) this corresponds to the continuous 
spectrum with threshold at ujq and density of states of the form l/{uj — ojq). 



Eq.(5J.) is also interesting from the following point of view. We expect that if the hadron at 
initial rapidity contains a small number of gluons, its wavefunction will have large projection mostly 
on the states in the Yang space, r]i ~ 0. Analogously if the hadron is very dense, it has a large 
projection onto the Yin space, 7^ ~ 0. Thus if we consider scattering of a dilute projectile on a dense 



target or vice versa, only one of the first two terms in the sum in eq.(5.1) contributes, while the 
third term is irrelevant. This corresponds to the situation described by the KLWMIJ or JIMWLK 
evolution. In this case one does not need to know the off diagonal matrix elements between the Yin 
and Yang eigenstates. The third term is important only if we are in the situation when both hadrons 
involved in the scattering are initially dilute. This is of course precisely the situation in which we 
expect the Pomeron loops to give important contribution to the evolution at high enough rapidity. 
Thus the last term in eq.(|5.lD is naturally associated with the contribution of Pomeron loops. 



This observation suggests an approximate way to include the Pomeron loops in the calculation 



of the scattering matrix via eq.(5T). Namely rather than searching for the full hamiltonian H^^^ 
one could solve for the eigenfunctions of }{J^^iwlk -^^ ^^^ yin space. Hi and for the eigenfunctions 
of ff^LWMij ^^ |.]^g Yang space, Gi. One can then use these two sets to approximate the third term 



in the sum of eq.(5.1). In spirit this is similar to approximations discussed in the recent literature 
||2l| , 23, 24, 12, ^] which strive to keep both the "Pomeron splittings" and "Pomeron mergings" in 
the evolution of the scattering amplitude. The advantage of the present suggestion is that it treats 
both processes symmetrically, and also keeps the essential nonlinearities in calculating each set of 
the eigenfunctions. It does not take into account any perturbative corrections to the eigenfunctions 
themselves, which will certainly be present in H^^^ . One may hope however that these corrections 
to the eigenfunctions Hi and Gi are indeed perturbatively small and negligible as long as the colliding 
particles at initial rapidity are perturbatively dilute. To pursue this route one would have to learn how 
to find the eigenfunctions and eigenvalues of the JIMWLK-KLWMIJ evolution, itself a formidable 
task. The potential payoff however seems to be worth investing some effort in this direction. 
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